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SYSTEMS  USING  COMPUTER  -  ALGEBRAIC  MANIPULATIONS 

by 

Rahul  Chattergy 
University  of  Hawaii 
Honolulu,  Hawaii 

Abstract 

J 

This  paper  discusses  the  reliability  of  operation  of  an  on-line  computer 
system  described  by  a  semi-Markov  process  model.  Analytical  solutions  are 
obtained  by  using  computer-aided  algebraic  manipulation  techniques.  This 
paper  demonstrates  that  the  difficulties  of  obtaining  analytical  solutions 
to  Markov  processes  by  standard  techniques  can  be  considerably  reduced  by 
the  application  of  algebraic  symbol  manipulation  languages.  To  the  author's 
knowledge,  the  results  of  the  reliability  analysis  are  also  new. 
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I.  INTRODUCTION 


Consider  an  on-line  computer  system  as  shown  in  Fig.  1.  The  users  of 
such  a  system  gain  access  via  the  terminals.  The  computational  needs  of  a 
user  are  satisfied  by  the  main  processor.  Communications  among  the  users, 
and  between  any  user  and  the  main  processor,  are  supervised  by  the  front-end 
processor.  The  main  purpose  of  a  front-end  processor  is  to  relieve  the  main 
processor  of  the  communications -management  duties.  In  such  a  system,  either 
the  front-end  processor  or  the  main  processor  or  both  may  fail  at  any  time. 

We  assume  that  if  the  front-end  processor  fails,  then  the  communications - 
management  job  is  taken  over  by  the  main  processor.  This  causes  an  increase 
in  response-time  as  observed  by  the  users  and  constitutes  a  degradation  of  the 
system-performance.  If  the  main  processor  fails,  then  evidently  the  users 
have  a  limited  computational  support  from  the  front-end  processor  but  they 
may  comnunicate  with  each  other  via  the  front-end  processor.  This  situation 
constitutes  a  more  serious  degradation  of  the  system-performance  than  the  one 
mentioned  above.  Finally,  if  both  processors  fail,  the  system  is  shut  down 
completely  for  repairs.  In  such  a  situation,  the  main  processor  has  priority 
on  repair  service  over  the  front-end  processor.  Note  that  the  system  is  not 
a  simple  special  case  of  the  repairman  problem  [2]  or  a  series  connection  of 
two  unreliable  subsystems  [3] . 

This  system  can  be  characterized  by  four  states  which  are  listed  in  Table 
1.  Transitions  between  pairs  of  states  occur  at  randomly  distributed  instants 
of  time.  Therefore,  it  is  possible  to  describe  the  system  by  a  semi-Markov 
The  reliability  of  operation  of  such  a  system  can  be  measured  by 
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various  probabilistic  measures  related  to  the  semi-Markov  process.  In  the 
next  section  we  describe  the  semi -Markov  process  in  more  detail. 


II.  SEMI -MARKOV  PROCESS  MODEL 

A  semi-Markov  process  is  described  by  a  set  of  states  and  probabilistic 
descriptions  of  the  temporal  transitions  between  pairs  of  states  [1] .  The 
on-line  computer  system  has  four  states,  each  state  uniquely  labelled  by  an 
integer  from  the  set  [1,...,4]  as  shown  in  Table  1.  Allowable  transitions 
between  pairs  of  states  are  shown  in  Fig.  2.  Let  t^  denote  the  time  spent 
by  the  process  in  state  i  before  a  transition  to  some  other  state  occurs. 

We  define  the  waiting-time  distribution  in  state  i  as 


wt(t)  =  ?[t±  <  t]  , 

and  the  corresponding  density  function  and  mean  by  w^ (t)  and  w^,  respectively. 

Let  [0 and  [0,F2]  respectively  denote  the  failure -free  operation  in¬ 
tervals  of  the  main  and  front-end  processors.  Following  argunents  in  [2]  and 
[3] ,  we  assume  the  failure-time  distributions  to  be 

P[F1  <  t]  =1  -  exp(-Ljt)  , 
and  P[F2  <_  t]  *  1  -  exp(-L2t)  , 

for  t>0  and  zeros  otherwise.  Let  [0,R^]  and  [0,R2]  denote  the  repair -time 
intervals  of  the  main  and  front-end  processors,  respectively.  We  assume  the 
repair-time  distributions  to  be 


i  t]  *  1  -  exp(-Gjt)  , 
P[R2  <_  t]  =  1  -  exp(-G2t)  , 


and 
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for  t>0  and  zeros  otherwise.  Note  that  since  the  main  processor  has  priority 
of  repair  service,  only  the  exponential  distribution  for  P[R2<t]  will  result 
in  a  Markov  process.  The  waiting-time  distribution  for  each  state  can  now 
be  evaluated  in  a  straightforward  manner.  As  an  example,  let  us  evaluate 
W^Ct).  We  define  the  following  events: 

E^:  Repair-time  of  the  main  processor  >  t  , 

E^:  Failure- time  of  the  front-end  processor  >  t  . 

Then 

1  -  W3(t)  =  P[Ej_  and  E2]  , 

=  ppy  P[E2]  , 

=  exp(-G^t)  exp(-L2t)  , 

=  exp(-(G1+L2)t)  ; 

and  hence 

W3(t)  =  1  -  exp(-(G1+L2)t)  . 

The  corresponding  density  function  is 

w3(t)  -  (G1+L2)exp(-(G1+L2)t)  ; 

and  the  mean  is 

w3  =  i/(g1+l2)  . 

A  list  of  the  four  waiting-time  distributions,  their  density  functions  and 
means  are  given  in  Table  2. 
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Let  p^j (t)  denote  the  conditional  probability  density  function  of  a 
transition  to  state  j  in  [t,t+A]  given  that  the  process  entered  state  i  at 

•ir 

i. 

time  zero  and  the  next  transition  from  state  i  occurs  in  [t,t+A]  for  suf¬ 
ficiently  small  A>0.  Then  the  probability  density  function  of  a  transition 

from  state  i  to  state  j  after  waiting  t  units  of  time  in  state  i  is  given  by 

+*  -ni 

• 

i 

cij  W  =  pij(t)wi(t)  • 

-4 

The  core  matrix  of  a  semi-Markov  process  is  defined  as  C(t)=[c— (t)] ,  and  it 

provides  a  complete  probabilistic  description  of  the  process  [1] - 

Let  us  consider  a  sample  calculation  for  the  conditional  transition 

-«r 

probability  density  function  from  state  3  to  state  2.  We  define  the  follow- 

.  •"  < 

! 

ing  events: 

E^:  Front-end  processor  fails  in  [t,t+A]  , 

V; 

9 

E^:  Main  processor  not  repaired  in  [0 , t]  , 

E,-:  Transition  occurs  out  of  state  3  in  [t,t+A]  . 

V 

Then 

i 

» 

P32OO  =  pKE3  and  E4)/E5]  , 

• 

-  _■< 

•  ■  i 

«  P[ (e3  and  E4)  and  E5J/P[E5]  , 

1 

=  P[E3  and  E4]/P[E5]  , 

”  ■. 

(AL2exp  (-d^t) )  exp  (-G^t) 

A  ( ( ]  +L  2 )  cx  p  ( -  ( c.  ^  +  L  2 )  t ) 

'1 

=  l2/(g1+l2)  . 

9 
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The  core  matrix  of  this  process  is  shown  in  Table  3.  The  single-step  trans- 

00 

it ion  probability  matrix  of  the  imbedded  Markov  chain  is  given  by  P=  /  P(t)dt. 

0 

The  P  matrix  for  this  process  is  shown  in  Table  4. 

Let  e^j (t) A  denote  the  probability  that  the  process  will  enter  state  j 
in  [t,t+A]  given  that  it  was  in  state  i  at  time  zero  and  let  E(t)=[eij (t)] . 
Then 

E(s)  =  [I-C(s)]  1 

where  I  is  the  identity  matrix  and  E(s)  and  C(s)  are  Laplace  transforms  of  E(t) 

and  C(t)  respectively  (see  {1]).  Define 

E  =  lim[sE(s)]  . 
s+O 

For  a  monodesmic  process  such  as  the  one  described  here,  the  rows  of  E  are 
identical  [1] .  Let  ej  denote  the  j  element  of  any  row  of  E.  Then  the  limit¬ 
ing  interval  transition  probability  for  state  j ,  denoted  by  in ,  is  given  by 

h.  =  e.w.  . 

J  J  J 

Suppose  the  process  has  been  operating  unobserved  for  a  long  period  of  time. 
Then  in  is  the  probability  of  the  event  that  the  process  will  be  in  state  j 
when  observed  next.  The  state  occupancy  statistics  can  also  be  obtained  from 
E(s).  The  details  of  this  and  mean  first  passage  time  computations  can  be 
found  in  [1] . 

The  next  section  shows  how  analytic  expressions  can  be  obtained  for  In , 
state  occupancy  statistics,  and  mean  first-passage  times  by  using  computer- 
aided  algebraic  manipulation  techniques.  This  material  was  generated  inter¬ 
actively,  by  using  the  symbol  manipulation  language  called  MACSYMA  available 


Table  1 


STATE  DESCRIPTION  TABLE 


STATE  LABEL 

Both  processors  in  operation  1 

Both  processors  down  2 

Only  main  processor  down  3 

Only  front-end  processor  down  4 


Table  2 

WAITING-TIME  DISTRIBUTIONS 


STATE 

DISTRIBUTION 

DENSITY 

MEAN 

1 

1  -  expOCLj+L^t) 

(L1+L2) exp ( - (L1+L2) t) 

i/(l1+l2) 

2 

1  -  exp(-G^t) 

G-jexpOG-^t) 

1/Gj 

3 

1  -  exp(-  (G^+L2)t) 

(G1+L2) exp ( - (Gx+L2) t) 

i/  (g1+l2) 

4 

1  -  exp(-  (L1+G2)t) 

(L1+G2)exp(-(L1+C2)t) 

l/C^+G^ 

E 


h 


Table  3 


GORE  MATRIX 

0  0  L^exp(-  (Lj+L2)t) 

0  0  0 
G1exp(-(G1+L2)t)  L2exp(-(G1+L2)t)  0 

G^exp  (-  (L,  +G-)  t)  L1expC*(L,+G-)t) 


L2exp(-  (L1+L2)  t) 
G^expC-G^t) 

0 


0 


0 
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Table  4 

TRANSITION  PROBABILITY  MATRIX  OF  THE  IMBEDDED  MARKOV  CHAIN 

o  o  l1/cl1+l2)  l2/ (l1+l2d 

0  0  0  1 
g1/(g1+l2)  l2/(g1+l2)  0  0 

G2/ai+G2)  L1/(L1+G2:)  0  0 


at  the  MIT  mathematics  laboratory,  accessed  via  the  ARPA  computer -comnum ca¬ 
tions  network. 


III.  INTERACTIVE  ANALYSIS  USING  MACSYMA 


Comments  imbedded  between  /*  and  */  are  not  part  of  the  interactive  session, 
but  have  been  added  to  explain  the  process. 


/*  TIME: TRUE  WILL  PRINT  CPU  TIME  USED  IN  EACH  STEP:  */ 
/*  A,B,C  ARE  CONSTANTS  USED  TO  SIMPLIFY  DATA  INPUT:  */ 

(Cl)  TIME: TRUE  $ 

TIME*  8  MSEC. 

(C2)  A: LI/ (L1+L2)  $ 

TIME=  35  MSEC. 


(C3)  B : G 1 / (G1+L2)  $ 

TIME*  28  MSEC. 

(C4)  C:G2/ (L1+G2)  $ 

TIME*  25  MSEC. 

/*  WAIT: ROW-VECTOR  OF  MEAN  WAITING  TIMES.  ELEMENT  I  */ 


/*  IS  THE  MEAN  WAITING  TIME  AT  STATE  I:  V 

(C5)  WAIT:MATRIX([A/L1 ,1/G1 ,B/G1 ,C/G2] ) ; 

TIME*  60  MSEC. 

[Ill  1  ] 

CDS)  [ . - . -  ] 


[  L2  +  LI  G1  L2  +  G1  LI  +  G2  ] 

/*  TCORE  IS  THE  CORE  MATRIX  OF  THE  SEMI-MARKOV  PROCESS  */ 

/*  TCORE  IS  PRINTED  ROW  BY  ROW,  EACH  ROW  ENCLOSED  IN  [  ] .  V 
/*  %E  DENOTES  THE  EXPONENTIAL  "e" .  “  DENOTES  EXPONENTIATION:  V 
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(C6) 


TCORE: MATRIX 

([0,0,L1*%E"((-L1/A)*T),L2*%Ea((-L1/A)*T)], 

[0,0,0,G1*%E~(-G1*T)J, 

[Gl *%EA ( (-G1/ B) *T)  , L2*%E~ ( (-G1/B) *T)  ,0,0]  , 
[G2*%E"(C-G2/C)*T),L1*%E"((-G2/C)*T),0,0]); 


TIME=  468  MSEC. 


(-  L2  -  LI)  T  (-  L2  -  LI)  T 

(D6)  MATRIX([0,  0,  LI  %E  ,  L2  %E  ], 

-  G1  T  (-  L2  -  Gl)  T  (-  L2  -  Gl)  T 

[0,  0,  0,  Gl  %E  ],  [Gl  %E  ,  L2  %E  ,  0, 

(-  LI  -  G2)  T  (-  LI  -  G2)  T 

0],  [G2  %E  ,  LI  %E  ,  0,  0]) 

/*  THE  NESTED  DO  LOOPS  COMPUTE  THE  LAPLACE  TRANSFORM  OF  V 
/*  TCORE  AND  STORE  IT  IN  A  TEMPORARY  ARRAY  TEMPO :  */ 

(C7) 

FOR  1:1  THRU  4  DO 

FOR  J: 1  THRU  4  DO 

TEMPO [ I, J]  -.LAPLACE [TCORE [ I, J]  ,T,S)  ; 

LAPLAC  FASL  DSK  MACSYM  BEING  LOADED 
LOADING  DONE 
TIME*  1195  MSEC. 

(D7)  DONE 

/*  GENMATRIX  OPERATION  IS  USED  TO  CONVERT  WRAY  TEMPO  */. 


/*  INTO  MATRIX  SCORE=LAPLACETRANSFORM  (  TCORE)  V 

/*  SCORE  MATRIX  SHOWN  BELOW  :  */ 


(C8)  SCORE: GENMATRIX (TEMPO, 4,4); 


GENMAT  FASL  DSK  MAXOUT  BEING  LOADED 
LOADING  DONE 
TIME*  98  MSEC. 

[ 

[  0  0 

[ 

[ 

r 

[  0  0 

f 

[ 

[  Gl  l: 

[  - . - . —  - 

[  s  +  r 2  +  Gl  S  ♦  L2  +  Gl 


LI  L2  ] 

. ] 

S  +  L2  +  LI  S  +  L2  +  LI  ] 

] 

Gl  ] 

0  .  ] 

S  ♦  Gl  ] 

] 
] 

0  0  ] 

] 
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t 

[  G2  LI 


[  S  +  LI  +  G2  S  +  LI  +  G2 


] 

] 

0  ] 
] 


(C9) 

TIME=  8  MSEC. 


KILL (TEMPO) ; 


(D9)  DONE 

/*  SENTRY=INVERSE  OF  (  IDENTITY  -  SCORE  ).  ~~  DENOTES  NONCOMMUTATIVE  */ 

/*  EXPONENTIATION.  INVERSE  OF  MATRIX  =  MATRIX  :  V 

(CIO) 


SENTRY:  (IDENT(4) -SCORE)''"-!  $ 


TIME=  82298  MSEC. 

/*  TEMPO=  S  *  FIRST  ROW  OF  SENTRY.  RATSIMP  SIMPLIFIES  EXPRESSIONS:  */ 

(Cll)  TEMPOS RATSIMP (S*ROW (SENTRY ,1))  $ 

TIME=  167759  MSEC. 

/*  LMT=  LIMIT  CF  TEMPO  AS  S—>0:  V 

(Cl 2)  LMT: LIMIT(TEMP0,S,0)  $ 

LIMIT  FASL  DSK  MACSYM  BEING  LOADED 
LOADING  DONE 
TIME=  2873  MSEC. 

/*  UiTDS2=  LIMITING  INTERVAL  TRANSITION  PROBABILITIES  */ 

/*  OF  THE  SEMI-MARKOV  PROCESS:  */ 

(Cl 3)  LMTDST: RATSIMP (LMT* WAIT)  $ 

TIME=  28120  MSEC. 

/*  CHECK  ON  LMTDST  FOR  L1=0  AND  L2=0.  EVALUATE  LMTDST  FOR  L2=0:  V 

(CI4)  LMTDST, L2=0  $ 

TIME=  683  MSEC 

(Cl  5)  RATSIMP (D14) ; 


TIME=  997  MSEC. 
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[  G1  LI  J 

(D15)  [  -  0  -  0  ] 

[  LI  +  Gl  LI  ♦  G1  ] 

/*  EVALUATE  LMTDST  FOE  L1=0:  V 

(Cl  6)  LMTDST, L1=0  $ 

TIME=  646  MSEC. 


(Cl  7) 


RATS IMP (D 16) ; 


TIME=  4652  MSEC. 


(D17) 


[  G2  L2  ] 

[  -  0  0  - ] 

[  L2  +  G2  L2  +  G2  ] 


/*  LIMITING  INTERVAL  TRANSITION  PROBABILITIES  IN  */ 
/*  GENERAL  FORM  SHOWN  BELOW  :  V 


(CIS) 


LMTDST; 


TIME=  2  MSEC. 

2  2 

(D18)  MATRIX ([(Gl  G2  L2  +  G1  G2)/((L1  +  Gl)  L2 

2  2  2 

+  (LI  +  (G~  2  Gl)  LI  +  Gl  G2  +  Gl  )  L2  +  Gl  G2  LI  +  Gl  G2), 

2  2  2 
(LI  L2  +  (LI  +  (G2  +  Gl)  LI)  L2)/((L1  +  Gl)  L2 

2  2  2 

+  (LI  +  (G 2  +  2  Gl)  LI  +  Gl  G2  +  Gl  )  L2  +  Gl  G2  LI  +  Gl  G2) , 

2  2  2 
Gl  G2  LI/ ((LI  +  Gl)  L2  +  (LI  +  (G2  +  2  Gl)  LI  +  Gl  G2  +  Gl  )  L2 

2  2  2 
+  Gl  G2  LI  +  Gl  G2) ,  (Gl  L2  +  (Gl  LI  +  Gl  )  L2) 


2  2  2 
/((LI  +  Gl)  L2  +  (LI  +  (G2  +  2  Gl)  LI  +  Gl  G2  +  Gl  )  L2  +  Gl  G2  LI 


2 

+  Gl  G2)] ) 

/*  COMPUTE  STATE  OCCUPANCY  STATISTICS.  T0CUP[l,2]  V 
/*  =  AVRG.  NS.  OF  TIMES  STATE  2  IS  VISITED  IN  [0}T]  */ 
/*  STARTING  IN  STATE  1  AT  TIME  ZERO.  ILT  =  INVERSE  */ 
/*  LAPLACE  TRANSFORM:  V 
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(C19)  S0CUP[1 ,2]  =  RATSIMP (SENTRY [1, 2] /S)  $ 

TIME=  5426  MSEC. 

(C20)  T0CUP[1 ,?]  =  ILT(S0CUP[1 , 2] ,S,T)  $ 

/*  INFORMATION  REQUESTED  BY  THE  ILT  ROUTINE:  V 

2  2  2 
IS  LI  +  (2G2  +  2G 1 ) LI  +  G2  -  2G1G2  +  G1 

POSITIVE,  NEGATIVE,  OR  ZERO? 

/*  ANSWER  ENTERED  FROM  THE  TERMINAL:  V 

POSITIVE; 

TIME=  162735  MSEC. 

(C21)  T0CUP[1 ,2] ,L1=1 ,L2=3,G1=10,G2=30; 

TIME=  5377  MSEC. 

(D21)  -  21  T  114998  SINH(SQRT(174)  T)  9515  COSH (SQRT (1741  T) 

%E  ( - ♦ - ) 

134657  SQRT(174)  134657 

-  6  T 

65  T  3  %E  841 


178  68  31684 

/*  NOTE  THAT  TOCUP[l,2]  HAS  A  LINEAR  TERM  IN  T  V 
/*  FOR  LARGE  I  r'HTS  TERM  WILL  BE  DOMINANT  :  V 

(C22)  PART (TOC UP [1,2] ,2); 

TIME*  3  MSEC 

2  2  2 

(D22)  (G1  LI  L2  +  (G1  LI  +  (G1  G2  +  Gl  )  LI)  L2)  T 

2  2  2 

/((LI  +  Gl)  L2  +  (LI  +  (G2  +  ?G1)  LI  +  G1  G2  +  G1  )  L2 

2 

+  Gl  G2  LI  +  Gl  G2) 


IV.  DISCUSSION  OF  THE  RESULTS 


Two  important  parameters  for  estimating  the  reliability  of  operation  of 
this  system  are  h^  and  l^,  respectively  the  probabilities  of  being  fully 
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2 

operational  and  completely  shutdown  in  the  steady  state.  L?t  D- (k2+G^G2) 
+(L^+(G2+2G1)L1+G1CG1+G2))L2.  Then 

hj  =  g1g2cg1+l2)/d,  (1) 

and  h2  *  L^L2(L2+L2+G^+G2VD*  (2) 

Equations  (1)  and  (2)  show  the  dependence  of  hj  and  h2  on  all  the  parameters 
of  the  system.  Let  us  investigate  the  effects  of  and  L2  on  hj.  Note  that 
L=0  means  the  corresponding  processor  never  fails.  As  L  goes  to  infinity, 
the  corresponding  processor  fails  with  probability  one  in  any  time  interval 
of  positive  length.  If  Lj=0,  then 

hj  «  G2/(L2+G2),  (3) 

-  1/(1+(L2/G2), 

1/ (1+ (Average  repair -time/ Average  failure-time)  ; 

and  if  L2*0,  then 

hx  =  Gj/CLj+G^,  (4) 

=  1/ (1+(L1/G1))  , 

=  1/(1+ (Average  repair-time/ Average  failure-time). 

For  either  L^*0  or  L2=0,  the  process  reduces  to  a  simple  two-state  semi-Markov 
process  and  h^  can  be  easily  computed  to  verify  the  above  results.  These 
results  were  used  as  checks  during  the  interactive  session  with  MACSYMA. 

As  L2  increases  without  bound,  (approximate)  and 

h,  =  G,/L.,  (approximate),  (5) 

X  L  i 

*  (Average  failure-time/ Average  repair -t  He)  for  front-end  processor; 


h2  =  L1/G1  (approximate). 


(6) 


=  (Average  reapir-time/Average  failure-time)  for  main  processor. 

2 

As  increases  without  bound,  D=LjL2  (approximate)  and 

2 

hj  =  G1G2(G2+L2)L1L2  (approximate),  (7) 

h2  *  1  (approximate).  (8) 

Let  us  now  consider  the  effects  of  G^  and  G2  on  h^ .  Note  that  G=0  means 
that  the  corresponding  processor  is  never  repaired,  and  as  G  goes  to  infinity, 
the  corresponding  processor  is  repaired  with  probability  one  in  any  time 
interval  of  positive  length.  As  G^  increases  without  bound,  I>G^ (G2+L2) 
(approximate)  and 

hl  =  G2^G2+L2^  (aPProxiraate^  •  (9) 

Note  that  equations  (3)  and  (9)  give  identical  values  of  h^  since  no  failure 
of  main  processor  is  equivalent  to  failures  followed  by  instantaneous  repairs. 
As  G2  increases  without  bound,  D=G^G2(G^+L^+L2)+G2L1L2  and 

hi  =  G^L^)  (approximate).  (10) 

Equations  (4)  and  (10)  give  identical  values  of  lv^  for  reasons  explained  above. 

Next  let  us  consider  a  state -occupancy  statistic  of  relevance  to  the 
reliability  of  this  system.  Let  N^CT)  denote  the  average  number  of  times  the 
process  visits  state  2  in  the  tune  interval  [0,T],  starting  in  state  1  at  time 
zero.  Then  for  large  values  of  T, 


Ni2(T)  =  G^T  (approximate)  . 


(ID 


Equation  (11)  is  derived  by  observing  the  results  obtained  interactively  with 
MACSYMA.  The  exact  expression  of  N^CO  contain  exponentially  decaying  terms 
of  T  and  the  linear  term,  shown  in  equation  (11).  It  is  interesting  to 

2 

consider  the  effect  of  on  N^CH.  As  increases  without  bound,  D=G^ (G^+L^) 
(approximate),  l^L^/CG^G 2+1.2))  (approximate)  and  consequently  G^h2=LjL2/ (G2+L2) 
Thus 

Ni2(T)  =  L^T/^)  (approximate).  (12) 

We  observe  that  L-^/ (G2+L2)  is  a  lower  bound  on  the  slope  of  N^CD  as 
goes  to  infinity.  Note  that  even  for  very  large  values  of  G^,  the  average 
number  of  times  a  complete  system-breakdown  occurs  in  [0,T]  cannot  be  made 
less  than  that  given  by  equation  (12). 

Next,  let  us  consider  the  average  first-passage  times  between  pairs  of 
states.  Let  T  denote  the  average  first-passage  time  from  state  i  to  state 
j ,  Then  we  define  as  the  recovery- time  of  the  system  from  complete  break¬ 
down  in  state  2  to  fully  operational  in  state  1.  The  unknown  quantities  T^ 
satisfy  a  set  of  linear  algebraic  equations  which  can  be  easily  solved  by 
MACSYMA.  Solving  these  equations  we  have 


T2i  =  (L^)  (1/Gj^  +  i/(l1+g2))/g2  . 


As  G.  increases  without  bound,  recovery-time  becomes 


Recovery- time  *  Average  repair-time  for  front-end  processor. 


As  increases  without  bound,  recovery- time  becomes 

T21  =  !/Gr 

Recovery- time  =  Average  repair-time  for  main  processor. 

These  two  values  of  T^  are  intuitively  obvious.  If  the  main  processor  does 
not  fail  too  soon,  i.e.,  is  close  to  zero,  then 

Tn  =  (1/G1)  ♦  (1/G2)  ,  (15) 

i.e..  Recovery- time  =  Average  repair-time  for  main  processor 

+  Average  repair-time  for  front-end  processor  . 

There  are  other  reliability -related  expressions  which  can  be  studied  in  a 
similar  manner.  We  believe  that  the  above  results  illustrate  the  usefulness 
of  computer-aided  algebraic  manipulation  techniques  for  analyzing  semi -Markov 
processes. 


V.  CONCLUSION 


The  main  advantage  of  using  computer-aided  algebraic  manipulation  stems 
from  the  fact  that  the  core  matrix  can  be  a  function  of  system  parameters. 

Closed -form  analytical  solutions  can  be  obtained  as  functions  of  these  system 
parameters.  Numerical  analysis  will  either  require  repetition  of  the  same 
computations  for  different  parameter  values  or  some  sort  of  sensitivity  analysis 
over  a  limited  range  of  parameter  values.  Also  there  is  no  need  to  invent 


ingenious  ad-hoo  methods  for  individual  problems.  Using  computer  algebra,  the 
same  standard  approach  can  be  used  to  analyze  all  such  processes. 

The  size  of  the  system,  as  measured  by  the  number  of  states,  that  can  be 
handled  in  this  manner,  obviously  depends  on  the  speed  and  memory-size  of 
the  computer  being  used,  It  has  been  observed  that  the  CRJ  time  required  to 
invert  a  matrix  to  compute  E(s)  is  usually  much  less  than  the  time  required 
for  simplifying  some  results  using  RATSIMP  or  the  time  required  for  computing 
the  inverse  Laplace  transform  to  obtain  (T) .  Normally  only  a  few  expressions 
of  interest  need  to  be  computed  for  any  system.  Hence,  as  long  as  E(s)  can  be 
computed  in  a  reasonable  period  of  time,  the  corresponding  system  can  be 
analyzed  by  this  approach. 
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